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Abstract 

We study the steady streaming between two infinitely long circular cylinders produced by small- 
amplitude transverse vibrations of the inner cylinder about the axis of the outer cylinder. The 
Vishik-Lyusternik method is employed to construct an asymptotic expansion of the solution of the 
Navier-Stokes equations in the limit of high-frequency vibrations for Reynolds numbers of order of 
unity. The effect of the Stokes drift of fluid particles is also studied. It is shown that it is nonzero 
not only within the boundary layers but also in higher order terms of the expansion of the averaged 
outer flow. 
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1. Introduction 

In this paper, we study two-dimensional oscillating flows of a viscous incompressible fluid between 
two infinitely long circular cylinders. The outer cylinder is fixed, and the inner cylinder performs 
small-amplitude harmonic oscillations about the axis of the outer cylinder. It is well-known that 
high-frequency oscillations of the boundary of a domain occupied by a viscous fluid not only produce 
an oscillating flow, but can also lead to appearance of a steady flow, which is usually called the 
steady streaming (see, e.g., [H). Here we are interested in this steady part of the flow. The basic 
parameters in our study are the inverse Strouhal number a and the dimensionless viscosity v (the 
inverse Reynolds number), defined as 

V* 1 v* 

a = ^ V= BTe = vfa (L1) 

where Vq is the amplitude of the velocity of the oscillating cylinder, a is its radius, uj is the angular 
frequency of the oscillations, and v* is the kinematic viscosity of the fluid. Parameter a measures 
the ratio of the amplitude of the displacement of the oscillating cylinder to its radius and is assumed 
to be small: a <C 1. One more dimensionless parameter which is widely used in literature is the 
'streaming Reynolds number' Re s = Vq 2 /lou* = a/u, which plays the role of the Reynolds number 
for the steady part of the flow. Note that for small a, Re s = 0(1) corresponds to v = 0(a), i.e. 
to high Reynolds numbers (Re = Ijv S> 1). 
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The steady streaming produced by an oscillating cylinder both in an unbounded fluid and in a 
cylindrical container had been studied by many researchers. A review of early works can be found 
in 0]. Regular perturbation analysis had been used to study both an oscillating cylinder in an 
unbounded fluid and the case of the flow between two cylinders in 0, 0] for the flow regimes with 
f>l (Re « 1). In [5j, the method of matched asymptotic expansions had been used to construct 
an asymptotic solution for an oscillating cylinder in an unbounded fluid for v = 0(1) (Re = 0(1)). 
In [a], the theoretical results of had been reconsidered, corrected with the Stokes drift, 

and compared with the experimental observations, which demonstrated a good agreement between 
the theoretical and experimental results. For high Reynolds numbers such that Re s = 0(1) or 
Re s 3> 1, the steady streaming induced by an oscillating cylinder in an unbounded fluid had been 
studied in 0, @, SH, flo| . For Re s = 0(1), the steady flow outside the Stokes layer is governed 
by the steady Navier-Stokes equations, and it had been shown that for Re s S> 1 there is a double 
boundary layer near the oscillating cylinder and that the steady streaming takes the form of a 
jet-like flow along the axis of oscillations. In 11, 12], the steady flow between two coaxial cylinders 



produced by small-amplitude transverse oscillations of the inner cylinder had been studied in the 
case of Re= = O(l) (and Re s 3> 1). The results of [121 ] are a good agreement with the experimental 



results of [13]. In [ll|, Il2j |. the case of small Re s corresponding to v = 0(1) had also been briefly 
discussed. 

The aim of this paper is to investigate steady streaming between two circular cylinders in the 
case of v = 0(1). As was mentioned above, this problem had been treated before in where 
an asymptotic solution of the Navier-Stokes equations had been constructed using the method of 
matched asymptotic expansions. Formulae given in can be used to write down a (composite) 
uniformly valid expansion for the averaged stream function. For an oscillating cylinder in an 
unbounded fluid, this had been explicitly done in [5]. However, the asymptotic solutions for both 
a single cylinder (presented in 0]) and for two cylinders (presented in [12]) are incomplete because 
of the absence of the 0(a 3 / 2 ) term associated with the averaged outer flow in the expansion of the 
averaged stream function (this will be discussed in more detail in Section 5). In the present paper, 
we construct a uniformly valid expansion of the averaged stream function in powers of a 1//2 up to 
terms of order 0(a 3 / 2 ). 

An interesting question which had not been discussed in 0, Qjl is the relation of the expansions 
obtained there to the original problem of a steady steaming produced by an oscillating cylinder as 
it is seen by an observer fixed in space. In [5j, the problem of a steady streaming produced by a 
cylinder which is placed in an oscillating flow had been considered. However, the steady component 
of the Eulerian velocity depends on the reference frame that is used to describe the flow. To clarify 
this point is one of the aims of this study. 



A similar problem arises in relation to the results of [12], where a conformal mapping that maps 
the gap between two eccentric cylinders onto the annulus between two cylinders with a common 
axis was employed. The subsequent analysis had been done using the transformed coordinates that 
depend on time, and the latter leads to a question whether a transformation back to the physical 
coordinates will change the steady component of the flow. The results of the present study give an 
answer to this question. 

Another question which is related to the previous two is the effect of the Stokes drift, which is 
understood here as a difference between the averaged Eulerian velocity field and the averaged 
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Lagrangian velocity (the velocity of fluid particles). The importance of Stokes drift is evident: (i) 
it is the Lagrangian velocity (rather than the Eulerian velocity) that is observed in experiments; 
(ii) it is the Lagrangian velocity that is invariant under the change of the frame of reference from 
the one fixed in the oscillating cylinder to the one fixed in space (and vice versa). The Stokes 
drift in various oscillating flows had been studied before by many authors (see, e.g., 14, H, 0, 0|)- 



However, there is still no clarity about when one can expect a nonzero effect of the Stokes drift in 
a flow produced by an oscillating cylinder. One of the aims of the present paper is to resolve this 
ambiguity. 

To construct the asymptotic expansion of the solution of the Navier-Stokes equations for small a, 
we employ the Vishik-Lyusternik method^ (see, e.g., 2, 17]) rather than the method of matched 



asymptotic expansions which is routinely used in fluid mechanics. In most cases that we know, at 
least first two terms in the uniformly valid asymptotic expansions produced by these two methods 
are the same. This does not mean that, given an expansion obtained by one of the two methods, 
it is easy to derive the same expansion using the other one. For example, although the asymptotic 
expansion constructed here using the Vishik-Lyusternik method can be obtained by the method 
of matched asymptotic expansions, this is not a straightforward procedure (it requires a certain 
transformation of the velocity field before the asymptotic procedure is started). In comparison 
with the method of matched asymptotic expansions, the Vishik-Lyusternik method involves more 
algebra in computing higher order terms, but it has two essential advantages: (i) it does not require 
the procedure of matching the inner and outer expansions and (ii) the boundary layer part of the 
expansion satisfies the condition of decay at infinity (in boundary layer variable) in all orders of 
the expansion, which is not the case in the method of matched asymptotic expansions where the 
boundary layer part usually does not decay and may even grow at infinity. The Vishik-Lyusternik 
method had been used to study viscous boundary layers at a fixed impermeable boundary by 
Chudov [18]. Recently, it has been applied to viscous boundary layers in high Reynolds number 
flows through a fixed domain with an inlet and an outlet [19|] and to viscous flows in a half-plane 



produced by tangential vibrations on its boundary [20(] 



In the present paper we compute first two nonzero terms in the asymptotic expansion of the steady 
velocity produced by an oscillating inner cylinder and the corresponding expansion for the stream 
function. The steady Eulerian velocity field is then corrected with the Stokes drift. The results for 
an oscillating cylinder in an unbounded fluid are obtained by passing to the limit R — )■ oo (where 
R is the ratio of the radius of the outer cylinder to the radius of the inner cylinder). 

The outline of the paper is as follows. In Section 2, we formulate the mathematical problem. 
In Section 3, we describe the method of constructing the asymptotic expansion and derive the 
equations and boundary conditions that are to be solved. These equations are solved and then 
corrected with the Stokes drift in Section 4. Section 5 contains the discussion of the results. 

2. Formulation of the problem 

We consider a two-dimensional flow of a viscous incompressible fluid between two circular cylinders 
with radii a and b (b > a) produced by small translational vibrations of an inner cylinder about 
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the axis of the outer cylinder which is fixed in space. Let x* = (x*,y*) be Cartesian coordinates on 
the plane and let Xq = (xo(i*), 0) be the position of the centre of the inner cylinder at time t*. We 
assume that is oscillating in t* with angular frequency uj and period T = 2tt/lv. Using 1/oj, 

a, Vq and pauoV^ as the characteristic scales for time, length, velocity and pressure (p is the fluid 
density), we introduce the dimensionless variables r = uf , x = ^x*, xq = ^Xq, v = y*v* and 



p = ■ The motion of the fluid is governed by the two-dimensional Navier-Stokes equations 

v r + e 2 [(v • V)v - i/V 2 v] = -Vp, V • v = 0. (2.1) 

Here e = y/a , a and v are the inverse Strouhal number and the inverse Reynolds number (the 
dimensionless viscosity), defined by (jl.ip . The velocity of the fluid satisfies the standard no-slip 
condition on the surfaces of the cylinders 

V ( X > T ) loutcr cylindcr= °> V ( X > T ) linner cylinder= ^ X o( T ) = ^o(r) e x . (2.2) 

Here dots denote differentiation with respect to r and xq(t) is a given function which prescribes 
the motion of the inner cylinder. In what follows we are interested in the asymptotic behaviour 
of periodic solution of Eqs. (|2.1|) . (|2.2|) in the high-frequency limit e — > 0. We assume that the 
amplitude of the oscillations of the cylinder is 0(e 2 ), i.e. xq(t) = e 2 f(r) for some 2-7r-periodic 
function / and / = 0(1). In what follows, we will consider f(r) given by 

/(r) = Re (Ce lT ) (2.3) 

where C is a complex constant having unit modulus (\C\ = 1). 
Boundary conditions (|2.2p take the form 

V ( X ' T ) louter cylinder = 0, v(x, t) | inner cylinder = / (r) e^. (2.4) 

The time-dependent boundary of the inner cylinder can be described in the parametric form by 
the equations 

x = cos 9 + e 2 /(r), y = sin § (2.5) 

where 9 £ [0, 2ir) is the parameter on the cylinder boundary. Now the boundary condition on the 
inner cylinder can be written as 

, =cose - +£ 2 /(T) = /'(r)e s . (2.6) 

i/— sin 

Using the assumption that e is small, we expand u and v in Taylor's series at point (x, y) 
(cos 9, sin 9) . This yields 



a; — cos O 



y— sin 



x—cos 



.+... = f'(r)e x (2.7) 



Note that each term on the left side of Eq. fj2.7f) is evaluated at the averaged position of the inner 
cylinder (where the axes of both cylinders coincide) . 
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In polar coordinates (r, 9) with origin at the axis of the outer cylinder, Eqs. (12. 1|) take the form 

,.2 / „ 2 



-Pr + e 



1 2 

r 



-u« r 1i6» H hf Vli 5 kVo 

V UV ( 2 u ^ 

— ItUr h z/ V f ^ + -^u e 



u 1 

u r -\ — h -f = o, 

r r 



(2.8) 
(2.9) 
(2.10) 



where u and v are the radial and azimuthal components of the velocity, subscripts V, V and '0' 
denote partial derivatives, and where V 2 = <9 2 + (l/r)d r + (l/r 2 )<9|. The boundary conditions at 
the outer cylinder are 



u 



r=R 



0, U 







(2.11) 



where R = b/a. The boundary condition (|2.7|) at the inner cylinder takes the form 

2 , r _ sin 6* 
+ e 2 / H v 

L r J r=l 

+ ... = /'(t)cos0, (2.12) 



r=l 



£ 4 /2 



+ 



2 sin 2 sin 2 sin cos 
L u u H k v 



r=l 



r=l 



Lv u 



r=l 



£ 4 /2 



+ 



L 2 u 



sin 2 2 sin I 



-Lu + 



2 sin cos 



-u 



-f(r) sm6. (2.13) 



r=l 



Here -L = cos 9 d r — (sin 9/r)de. 



3. Asymptotic expansion 



We seek a solution of (|2.8|) - (|2.13|) in the form 

u = u\r, 9, r, e) + e u a {£, 9, r, e) + e n b (??, 0, r, e) 
w = ^(r, 0, r, e) + v a (£, 9, r, e) + v b ( V , 9, r, e), 
p = p l (r, 9, r, e) + p a (£, 9, r, e) + p\r } , 9, r, e). 



(3.1) 
(3.2) 
(3.3) 



Here £ = (r — l)/e and = (R — r)/e are the boundary layer variables. Functions u l , v l , p % 
represent a regular expansion of the solution in power series in e (an outer solution), and u a , v a , 
p a and u b , v b , p b correspond to boundary layer corrections to this regular expansion. Superscripts 
'a' and 'b' correspond to the boundary layers at the inner and outer cylinders respectively. We 
assume that the boundary layer part of the expansion rapidly decays outside thin boundary layers, 
i.e. u a ,v a ,p a — > as £ — > oo and u b , v b ,p b — > as rj — > oo. 

3.1. Regular part of the expansion 

To derive the equations governing the regular part of the expansion, it is convenient to work 
with the Cartesian form ()2.ip of the Euler equations. Later we can rewrite the results in polar 
coordinates. Let 



= vj + e vi + e 2 v 2 + . . . , p* =p\ j + ep\ + e 2 p| + . . . 



(3.4) 



where v* = v\e x + v\e y and v k = v\ k e x + v l 2k e y (k = 0,1,2,...). On substituting (|3.4|) in 
()2.ip and collecting terms of equal powers of e, we find that the successive approximations v^, p' k 
(k = 0, 1, 2, . . .) satisfy the equations: 



d T v\ = -Vpl, V-v* fc = (3.5) 



for k = 0, 1 and 



d T v\ = -Vp\ + (- £(vj • VK_ 2 -z + ^V 2 v^ 2 ^ , V • v* fc = 0. (3.6) 
for k > 2. In what follows, we will use the following notation: for any 2-7r-periodic function F(t), 

F(t) =F + F(t), F=±-j F(r)dr, 

o 

i.e. F is the mean value of F(t) and, by definition, F(r) = F(t) — F is the oscillating part of F(t). 



Leading- order terms. Consider first Eqs. (|3.5p for k = 0. We seek a solution Vq which is periodic 
in t. Averaging the equation for Vg yields: Vpg = 0, which is the necessary condition for existence 
of periodic (in r) solutions for v . Without loss of generality, we put = 0, i.e. in the leading 
order the pressure in the outer solution is purely oscillatory with zero mean value. The general 
solution of Eqs. (|3.5[) can be written as Vq = Vq + vo where v = V</>o and (fio has zero mean value 
and is the solution of the boundary value problem 

V 2 0o = O, <p 0r \ r=1 = f'(r) cos 9, (f> 0r \ r=R =0. (3.7) 

The boundary conditions for <j)Q r at r = 1 and r = R will be justified later. 

On averaging the equations for v?> and the second of equations ()3.5|) . we obtain 

(vj • V)v^ + (v • V)v* = -Vpl + uV 2 %, V-v = 0. (3.8) 

Since v is irrotational, we can rewrite (|3.8p as 

(vj • V)v = -Vn + ^V 2 v , V • v = 0, (3.9) 



where IIo = + l^^o| 2 /2. Equations (|3.9p represent the time-independent Navier-Stokes equations 
that describe steady flows of a viscous incompressible fluid. Boundary conditions for v will be 
specified later. 

First-order terms. The solution of Eqs. (|3.5p for k = 1 has the form = v\+v\ where = V0i 
and (j>\ has zero mean value and is the solution of the boundary value problem 

V 2 0i = O, <p lr \ r=1 = ai (8,T), cf> lr \ r=R =bi(e,r). (3.10) 

Functions cl\{0,t) and bi(9,r) will be defined later. Manipulations similar to those employed in 
derivation of Eqs. (|3.9p lead to the following equations for v|: 

K-v)vi + (v i 1 -v)v^ ) = -vn 1 + ^v 2 vi, v-vi = o, (3.11) 



where Hi = p l 3 + (V0o • V0i). Boundary conditions for v\ will be specified later. 
Second-order terms. It will be proved later that 

= and v\ = 0. (3.12) 
Using these and the fact that both Vq and v\ are irrotational, it can be shown that 

d T v\ = -VQ 2 , V • v\ = 0, 

where Q2 = p\ + (^^°) _ it follows that v?, is irrotational, i.e. yr\ = V</>2, and (f>2 is the 

solution of the boundary value problem 

V 2 2 = O, </> 2r | r=1 =a 2 (0,r), ^2r| r=fi =6 2 (0,r), (3.13) 

where functions 0,2(0, r) and 62 t) will be defined later. 
The equations for v| can be written in the form 

= -Vn 2 + i/VV 2 , V-v| = 0, (3.14) 

where II2 = p\ + V^o • V</>2 + ^ v< ^ and where we have used the assumptions (|3,12p . Thus, the 
second order averaged velocity is described by the Stokes equations. Again, boundary conditions 
for v| will be specified later. 

Third-order terms. Separating the oscillating part of the equation for V3 and employing (|3,12p and 
the fact that Vq and v\ are irrotational, we find that 

5 T v| = -VQ 3 , 

where Q3 = p 3 + V0o • V0i — V(/>o ■ V0i. This equation and the continuity equation V • V3 = 
imply that v 3 = V(p 3 where 4> 3 is the solution of the boundary value problem 

V 2 ^3 = 0, 03r| r= i= a 3 (9,T), (t>3r\ r=R =b 3 (6,T). (3.15) 

Function a 3 (9,r) and b 3 (6,T) will be defined later. 

The equations for v| can be derived in the same manner as the equations for Vq, vj and v 2 . They 
are given by 

= - Vn 3 + ^V 2 v|, V-v| = 0, (3.16) 
where 1I3 = p| + V</>o • V^3 + V^i • V<p2- 

Thus, we have found that both the second and third order averaged velocities (v| and v|) in the 
outer flow are solutions of the Stokes problem with boundary conditions which will be determined 
later. 
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3.2. Boundary layers 

Boundary layer at the inner cylinder. We assume that 

u° = «g + e«f + ..., v a = v% + ev{ + p a = p a + ep\ + . . . . (3.17) 

Now we use our assumption that u b , v b and p b are nonzero only within a thin boundary layer 
near the outer cylinder and drop them from Eqs. ()3.1[) — f)3.3[) . Then we substitute the resulting 
equations, as well as (j3.4[) and (|3.17p . into Eqs. (|2.8|) - (|2.10p and take into account that u k , v k , 
p\ (k = 0,1,...) satisfy the equations (|3.5p . (|3.6p . After that, we make the change of variables 
r = l + e£, expand every function of in Taylor's series at e = and collect terms of equal powers 
in e. This yields the following equations: 

w2r-Ke*+Pfc* = *k, (3-i8) 

Pk = G a k , (3.19) 
u a H + v a ke = Hl (3.20) 

for k = 0, 1, . . . Here F a = 0, Gg = and Hg = 0; for k > 0, functions F fe a , G£ and H% are defined 
in term of v , . . . , v l k _ 1 , Ug, • • • , v k _ l and ttg, . . . , u k _ 1 . Explicit expressions for these functions are 
given in Appendix A (for k = 1,2, 3). 

Boundary layer at the outer cylinder. Let u b = Uq + eu b + . . ., v b = Vq + e v b + . . . and p b = 
p b Q + ep\ + ... The same procedure as before produces the following sequence of equations: 

v b kT - vv\ m + ^p b k9 = F b , (3.21) 

Pin = °l ( 3 - 22 ) 
-4 v + ^ b k e = H b k , (3.23) 

for k = 0, 1, . . . Here Fq = 0, Gq = and Hq = 0; explicit expressions for F k , G k and H k are given 
in Appendix A (for k = 1, 2, 3). 

We require that in all orders the boundary layer corrections to the outer solution rapidly decay 
outside boundary layers, i.e. (for each k = 0, 1, . . .) 

u%, v k , p k — > as £ — > oo and v b k , p\ — > as r/ — >■ oo. (3.24) 



3.5. Boundary conditions 

Now we take our expansions of the velocity in the outer flow and in the boundary layers, substi- 
tute them into f)2.11j) — (|2.13j) and collect terms of equal powers in e. This produces the following 
boundary conditions: 



4lr=i +«fc-iU=o = Qk, (3-25) 



4l r=1 +^L =0 = S k , (3.26) 



4lr=i?+4-ll„=0 = °> ( 3 - 27 ) 
.6 



4l r=jR +^|^o = 0. (3.28) 



for k = 0, 1, . . . Here 

Qo = f'if) cos 0, S = -f'{r) sin 9, 

functions Qk and Sf* for k > depend on Vq, . . . , ~v\_ 1 , Vq, . . . , v\_ x and Uq, . . . , and are given 
in Appendix A (for k = 1, 2, 3). Note that boundary conditions for <j)Q r at r = 1 and r = R in the 
boundary value problem (|3.7|) follow directly from (|3.25|) and (|3.27|) with A; = 0. 



4. Analysis of the asymptotic equations 

4-1. Leading order equations 

Outer flow. The solution of (|3.7p that describes the (leading order) oscillating outer flow is 

^ = -W^i{ r + T) CCB0 - (41) 

Inner cylinder. Consider now Eqs. (|3.18p ~ (|3.20p for k = 0. The condition of decay at infinity (in 
variable £) for p^ and Eq. f|3. 19j) have a consequence that = 0- Equation (|3. 18|) simplifies to the 
standard heat equation 

v a 0r = wfo (4.2) 
Boundary condition for v q at £ = follows from (|3.26|) (with A; = 0): 

^ol 5 =o= "«olr=i = -^— x Re{iCe lT ) sinO. (4.3) 

The solution of (|4,2p subject to the boundary conditions (|4,3p and ()3.24p is given by 

v a Q = -^^Re(iCe-< + *)sme, 7 = ^- (4-4) 

It follows from (|4.4p that Vq = 0. Thus, in the leading order the boundary layer at the inner 
cylinder is a purely oscillatory Stokes layer. This fact implies that the boundary condition for v 
at r = 1 (that is obtained by averaging the condition (|3.26|) ) is v l \ r=1 = 0. Similarly, averaging the 
condition (|3,25p yields u \ r=1 = 0. Thus, we have 

v£| r=1 =0. (4.5) 



The normal velocity Uq is determined from Eq. (j3.20p : 

oo 

«8&0,t) = J vfotf,6,T)d? = - -^—^ (^Ce-^A co S 9. (4.6) 

Here the constant of integration is chosen so as to guarantee that Uq(£,6,t) decays as £ — > oo. 
Uq|^ =0 gives us the boundary condition for the next approximation of the outer solution. Indeed, 
according to (|3.25[) for k = 1, we must have 

u\(r,6,r)\ r=0 = -u8(60,t)| c=o = J^;^ f^ 6 ") C ° s6 > (47) 

9 



This equation defines function a\{0,r) in (|3. l(Jj) . 

Outer cylinder. Consider now Eqs. (|3.21j) — f)3.23|) for k = 0. An analysis similar to what we did for 
the boundary layer at the inner cylinder results in the formula 

vl = --JL_ Re [iCe-^ +iT ) sin 9. (4.8) 
As before, the radial velocity Uq is determined from the incompressibility condition f|3.23j) : 

oo 

u o = -^J vM,0,T)d V ' = R{R 2_ 1) Re (~Ce-rt +i A cosfl. (4.9) 
v 

Again, the constant of integration is chosen so as to guarantee the decay of Mq as r\ — > oo. WqI^q 
gives us the boundary condition for the next approximation of the outer solution: 

U\ | r=1 = -U b Q \ v=0 = - RiR 2_ 1} Re {^ CelT ) COS6 - C 4 - 10 ) 

This equation defines function bi(9,r) in (|3.1Up . 

It follows from (|4.8p that Vq = 0, i.e. in the leading order the boundary layer at the outer cylinder 
is purely oscillatory. This, in turn, implies that the boundary condition for v l at r = R (obtained 
by averaging the condition (|3.28p ) is Vq\ t=r = 0. Similarly, Eq. ()3.27|) yields UQ\ r=R = 0. Hence, 

v*| r=i? =0. (4.11) 

Averaged outer flow. Equations (|3.9p and boundary conditions ()4.5p and (|4.1ip imply that Vq = 0, 
i.e there is no steady streaming in the leading order of the expansion. This justifies our earlier 
assumption about Vq. 



4-2. First order equations 

Oscillatory outer flow. Since now we know functions ai(9,r) and bi(6,r) (defined by Eqs. ()4.7p 
and (|4.10p ). we can solve problem (|3.1U|) . The solution is given by 

W ( {R + l)r + Re ( —eA cosO. (4.12) 



(i? 2 - l) 2 V r J V 7 

Inner cylinder. Consider Eqs. (l3~T8l) - (^20l) for k = 1. The fact that Gf = (Eq. (jA~4l in 
Appendix A) and the same arguments as before lead us to conclusion that pi = 0. Hence, Eq. 
(|3.18p reduces to 

< ~ = " cos 9 f'(T)vfc + uv^. (4.13) 
Averaging in r and integrating in variable £ twice, we find that 



i cos6< J f'(T)v$d£'. (4.14) 



v 1 -- 

is 
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Here the constants of integration are chosen so as to ensure that vf — > as ^ — > oo. It follows from 
the definition of averaging that / 1 (r)/2(r) = — f\{r) f' 2 {T) for any 27r-periodic functions /i and ji- 
Employing this property in (|4.14p . we obtain 

oo oo 

v\ = ~ cos0 J WW T d£ = co S I f(r)v a m d£ = -cos#70>^- (4.15) 

Here we used the fact that v$ satisfies Eq. (fPj) . Substitution of (f2T3j) and ([O} in (|4.15|) yields 

-5/V2^ (r^ltUFhA - ^(PU/o7,\\ em 9fl (4.16) 



cos(^/V2^) - sin(f/V2i/) sin 20. 



1 2y^ ^ 2 " 1 
The oscillatory part of i?" satisfies the equation 

"it - = ~ cos 6 f' v k + vv k ' 

the condition of decay at infinity and the boundary condition 



2*tf + * + 2) /iC, 



(i? 2 - l) 2 



7 



sin — cos 9 fv?:, 



0£ 



(4.17) 



(4.18) 



which follows from the oscillatory part of (|3.26|) and from Eq. (|4.12|) . Standard but tedious 
calculations result in 



cos 9 fv^ +w, w 



2R 2 



R 2 -l 



Re 



iC 



(R? + R + 2 1 £ 
\R(R 2 -l) 7 2," 



'7? + «T 



8111 ( 



(4.19) 



We do not give an explicit formula for fv^ as we do not use it in what follow^]. Both uf and u\ 
are computed using Eq. (|3.2Up . We have 



R 2 



R 2 -l 



sin^/V^) cos 29, 



d e (cos0M) +h, h 



2R 2 



R 2 -l 



Re 



C fR 3 + 3R + A 1 3£ 
R{R 2 -l) 2^~ ~2 



Outer cylinder. Similar analysis, applied to Eqs. (|3.2ip — (|3.23f) . yields 



2 sin l 



R{R 2 - 1 
2cos6> 
R 2 (R 2 -1) 



-Re 

i 

Re 



iC 



2R 3 + R 2 + 1 1 



R 2 -l 



7 2 



iC (4R 3 + 3R 2 + 1 1 3r? v 
"7 V B 2 ^ 27 + "2 1 ' 



(4.20) 
cos0. (4.21) 

(4.22) 
(4.23) 



Equations (|4T22l) and (j4T23j) imply that ^ = and u\ = 0, i.e. in contrast with the inner cylinder, 
there is no first- order steady boundary layer at the outer cylinder. 



fv^ as a function of r is purely oscillatory with double frequency of oscillation and therefore, it produces zero 
contribution to all quantities which will be of interest to us. 
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Averaged outer flow. On averaging boundary conditions (|3.27|) and (|3,28|) (for k = 1) and using 
the fact that tig = and v\ = 0, we find that 

*\l=R= 0- (4-24) 
Further, averaging boundary conditions (|3.25|) and (|3.26p and using (j4.6[) and (|4.15p . we obtain 

«llr=l= -^0^=0= °> V\\r=l= -Wllf=Q " cos 6 f( T H^=0 = °- 

These, together with (|4.24p . imply that Eqs. (|3.1ip should be solved with zero boundary conditions, 
which, in turn, leads to a conclusion that v\ = 0. This justifies our earlier assumption and means 
that there is no steady outer flow in the first order of the expansion. 



4-3. Second order equations 

Oscillatory outer flow. To find the oscillatory part of the second-order outer flow we need to solve 
problem (|3.13p for the velocity potential fa. Functions 0,2(6, t) and 02(6, r) which appear in (|3.13p 
are determined by boundary conditions (|3.25p and (|3.27p for k = 2. With the help of Eq. (|3.25p 
with k = and the continuity equations for vg and v , boundary condition (|3.25|) for k = 2 can 
be reduced to 

4l=i= -«il e =o +°e ( sin 6 Hl=i) = -«il e =o + cos 26 //'■ ( 4 - 25 ) 

The oscillatory part of (|4.25p gives us 0,2(6, r): 

o 2 (6,t) = -fi?| f=0 + cos 26 Jf. (4.26) 

The second term on the right side of (|4.26p will be ignored because it makes zero contribution to 
all quantities which we are interested in. The oscillatory part of (|3.27p yields 02(6, t): 

b 2 (6,r) = -u b 1 \ v=0 . (4.27) 

Substituting ()4.2ip in (|4.26p and (|4.23p in ()4.27p (and ignoring the second term on the right side 
of (|4.26|) as well as the first term in (|4.21|) ). we solve problem (|3.13|) . The solution is given by 

/ cos 6 ( +l?r+ ^-2rf» + 6tf-2fl + l\ _ 



(R 2 -l)(R-lf 



Inner cylinder. Consider now Eqs. (|3. 18p — (|3.2Q|) for k = 2. First we note that it follows from Eqs. 
(|4.2p and (|4.6p that Uq t = vuq^, which, in turn, implies that G% = (see Eq. (|A.5[> in Appendix 
A). Hence, Eq. (|3.19p reduces to = 0. This equation and the condition of decay at infinity 
imply that p\ = 0. Equation (|3. 18[) takes the form 



j 2t 



uv^ = F%. (4.29) 



Averaging yields the equation = —\F^ where F§(^,6) is obtained by averaging Eq. (|A.2p . 
The solution of this equation that satisfies the condition of decay at infinity can be written as 



00 00 

1 



F2(£",0) <%"<%'. (4.30) 



12 



Lengthy, but standard calculations result in the formula 
_ a Re~hin(29) 



A(R 2 - l) 2 



2R 3 e~ i - (5R 3 - R) £ (cos £ - sin£) + 4i? 3 cos £ + (lli? 3 + R - 4) sin £ 



(4.31) 



where £ = £/v2za The oscillatory part of can be obtained by separating the oscillatory part 
of Eq. (|4.29p and solving it. We will not do it as we are only interested in the steady part of the 
solution. Equation (|3.20[) for k = 2 is used to obtain u\'- 



^Rertcos 29 
A(R 2 - l) 2 



2R 3 e~t + (14.R 3 - 6i?) £sin£ + (12R 3 - 4) sin£ + (20i? 3 - 4) cos£ 



Outer cylinder. Similar calculations result in the following formulae for v% and u\: 



"2 
-b 



e v sin 20 r _f, -, 
2ft(ft2 iy2 l e 2fj (cos fj - sin fj) + 8 sin fj + 2 cos fj\ 



yJ2v e-v cos 29 r _= A 1 
u 2 = 2 R 2 (R 2 - l) 2 + 4r / smr ? + 8smr / + 12cosr /J > 



cos 29. 
(4.32) 

(4.33) 
(4.34) 



where fj = r)/y/2v. 

Averaged outer flow. The steady part of the second-order outer flow is determined from the Stokes 
equations (|3.14p . Averaging (|4.25p and using yields (|4.20p . we obtain 

4lr=l= °" ( 4 - 35 ) 

Boundary condition for v\ is obtained by averaging ()3.26p . It can be shown that it reduces to 

3 R* 



V 2\r=r 



2 (R 2 - 1)' 



■sin 29. 



Similarly, it can be shown that 



u 



2\ r =R 
4\ r =R 



n 1 77=0" 



- V 2\r)=0 z 



0, 



2R(R 2 -iy- 



sin 20. 



(4.36) 

(4.37) 
(4.38) 



It is convenient to introduce stream function ip\ such that u\ = and v\ = —ip\r- Then the 
Stokes equations (|3. 14[) and the boundary conditions (|4.35p - (|4.38p take the form 



v 4 vi = 0, n r =i= n r=R - "1 

- 3 RA • « fl 7t , _ 3 1 

2 ' lr=1 ~ 2 (R 2 - 1)2 Sm ^' ^2rl,= R - 2 fi(jR2 _ 1)2 



0, 



sin 2(9. 



The solution of (|4.39p is given by 



■n 



4(R 2 - l) 4 



Ci + ^ + C 3 r 2 + C 4 r 4 I sin 20 



(4.39) 



(4.40) 



where Ck (k = 1, . . . , 4) depend only on R and are given in Appendix A (Eq. fA25]) \ Equat ion 
(|4.40p is in agreement with earlier results of Duck and Smith (Eq. (3.33) in [12|) and Haddon and 
Riley (Eq. (3.1) in [111]). Formulae for u\ and v\ can be easily obtained from Eq. (|4.40p . 
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4-4- Third order equations 

Here we are interested only in the averaged outer flow. However, to find it, we need boundary 
conditions which come from the averaged boundary layers. 

Inner cylinder. Averaging Eqs. (|3.18p and (|3,19p for k = 3, we get 

1 „ 1 



v 



■3« = ~i&--*3 and ^ = G a 3 , (4.41) 



where Fg and G3 are obtained by averaging Eqs. ()A.3p and (|A.6p . We do not need explicit solutions 
of Eqs. ()4.4ip . all we need is VjjI^q. First, we integrate the second equation (|4.4ip . Then we insert 
P3 into the first equation ()4.4ip and integrate it twice. This yields 



00 00 00 00 00 

1 j J f - a . m m n . 1 



^ = --J J J G 3*(r , 0) di"' di" d?--J J F 3 a (e", 0) di" (4.42) 

Outer cylinder. Similarly, it can be shown that 

00 00 00 00 00 

4 = ™£ / / / &M",<>) H" dv" dr l'-\J J H(V", 0) H' drj. (4.43) 

V r)' rj" rj if 

Here F\ and G\ are obtained by averaging Eqs. (|A.12p and (]A.15p . 

Averaged outer flow. The steady part of the third-order outer flow is a solution of (|3.16p subject 
to appropriate boundary conditions. These boundary conditions are obtained by averaging (|3.25|) - 
(|3.28p for k = 3. Averaging Eqs. (|3.25p and (|3,27p yields the equations u\\ r=1 = — -uf |^ =0 an d 
=0 . Substituting (|4.32p and (|4.34p into these equations, we obtain 

«3lr=!= 2 (^_1)2 CQS2 ^ U ^=R= 2R 2 {R 2_ 1) 2 COs2g - ( 4 " 44 ) 

On averaging Eq. (|3.26p for k = 3 and using zeroth- and first-order boundary conditions, the 
boundary condition for v\ at r = 1 can be simplified to 



'3lr=r 



f2 f3 
,2 q£ „.a I „ A „3 a J 



^=0 -cosOf [y\ r \ r=1 +*>2$l$=oJ - cos d Y V ^^=°~ COS 9 ~Q V ^k=o- 
Substituting here Eq. (|4.42p and the explicit formulae for Vq, v\, Uq, u\, Vq, vf, etc., we find that 
/— , N , % (53i? 4 - 53i? 3 - 12i2 2 + 4fi - 4) 

Similarly, by averaging (|3.28[) and using (|4.43p . we obtain 

4U R =-^K 2 (R) S m2e, K 2 (R) = ^^R-f)^\)* ~ ' {AM) 
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On introducing stream function ip\ such that u\ = \^\q an d v\ = —ip\ r , Eqs. (|3.16|) and boundary 
conditions (|4.44|) - (|4.46|) can be written as 



V 4 ^ = 0, 



^i\ r=l =-V2^K 1 {R)shi2e, $Ur=R= V2vK 2 {R)sin2e. (4.47) 
The solution of this boundary value problem is given by 

* " WR-»iB>-W { Dl + ^ + Dp3 + D>r> ) Si ° 26 <448) 

where (k = 1, ... ,4) depend on R only and are given by Eq. (|A.26|) in Appendix A. Explicit 
expressions for u\ and v\ can now be easily obtained from (|4.48[) . 

4-5. Stokes drift 

So far we have discussed the Eulerian velocity. However, the velocity observed in experiments is 
the velocity of fluid particles, i.e. the Largangian velocity. It is well-known that in oscillatory 
flows the observed averaged Lagrangian velocity differs from the averaged Eulerian velocity, and 
the difference between these two is known as the Stokes drift. Below we discuss the effect of the 
Stokes drift on the averaged Largangian velocity. 

The motion of fluid particles is governed by the ordinary differential equation 

^ = 6 2 v(x,r,e), (4.49) 

The velocity field v(x, r, e) is the solution of the Navier-Stokes equations, which is 27r-periodic in 
r and has a nonzero average. We have already computed first three terms in the uniformly valid 
asymptotic expansion of v(x, r, e). Now we are interested in constructing an asymptotic expansion 
of the solution of (|4.49[) for small e. We introduce the slow time t = e 2 T, assume that x = x(r, t, e) 
and substitute this in f|4.49j) . This gives us the equation 

x T = e 2 (v(x,r, e) - x t ) . 

It is shown in Appendix B that the solution of this equation can be presented in the form x(xq, t) = 
x(xo,T,t,e) + 5t(xo,t,e) where x(xo,r, t, e) represents purely oscillatory part of the motion of the 
fluid particle whose averaged position at t = was xo and x(xo, t, e), is the solution of the equation 

x t = v(x,e) + v s (x,e). 

this equation describes slow motion of this particle due to the steady part of the Eulerian velocity 
field, v(x, e), and the Stokes drift velocity, v s (x, e). If we denote the Lagrangian velocity of fluid 
particles by superscript L and the Eulerian velocity by superscript E, then our results can be 
summarized as follows. Our asymptotic expansion for the averaged Eulerian velocity has the form 



u E = e 2 [4 + ui] + 0(e 3 ), v E = evf + e 2 v\ + v% + v\ + 0(e 3 ). (4.50) 
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It is shown in Appendix B that the Lagrangian velocity of fluid particles is given by 



it 



v\ + + v\ + v s 2 



+ 0(e 3 ) 



(4.51) 



where the Stokes drift velocity of the fluid particles v 2 is given by Eq. (fR33j) in Appendix B. 
Comparing (|4.50p with (|4.51j) . we observe that the Stokes drift eliminates (i) uf from the first 
equation (|4.50p and (ii) the O(e) term from the second equation (|4.50p . It also results in the 
additional 0(e 2 ) term v 2 in the expansion of the azimuthal velocity. Thus, the 0(e) steady boundary 
layer at the inner cylinder disappears when we take account of the Stokes drift. This is a consequence 
of the fact that the steady Lagrangian velocity rather than the steady Eulerian velocity are invariant 
to the change of reference frame. More precisely, the steady Eulerian velocity is not invariant in 
the following sense: it can be shown that in the reference frame fixed in the inner cylinder, there 
would be no 0(e) steady boundary layer at the inner cylinder and an O(e) steady boundary layer 
would appear near the outer cylinder (which is oscillating in this reference frame). It can be also 
shown that the steady Lagrangian velocity is the same both in the oscillating reference frame and 
in the fixed one. We note here that if one uses the method of matched asymptotic expansions, then 
although it is possible, it is not easy to detect the existence of an 0(e) steady boundary layer near 
the inner or outer cylinders. Fortunately, the presence or absence of the O(e) steady boundary 
layer do not affect the outer flow: the O(e) steady outer flow is zero in both cases. 

Further calculations with the help of the known formulae for v 2 , v 2 and v 2 show that v can be 
written as 

(4.52) 



v l 2 + (v a 2 ) L + (4) L + 0(<?) 



where 



R A sin 29 



F'(0, (4) 



sin 29 



4(^-1)2 yi> AR(R 2 -1) 2 

where F(z) = e~ z (3e~ z + 8cosz + 8sin2:), £ = £/y2i> and r) = r\j\f2~u. 



(4.53) 



4-6. Asymptotic expansion for stream function 

To rewrite our asymptotic expansion of terms of the averaged stream function, we first observe 
that 



+ e z 



where tp k is such that u\ = and v\ 



-ip l kr for k 



+ e° 



V<3 + i>2 + M 



+ 0(e 4 



0, 1, . . . and where ip k , ip b k are defined as 



, -j v a k (£,e,T)dg, ^ 



v b k (ri',9,T)d V ' 



for k = 0, 1, . . . Similarly, we have ip = e 2 Uf)\ + ipf] + e 3 [V'l + V'f + ^2] + and 

^ = € 2 fy + £ 3 r ( ^)L + ^a )L + ^L] + 0(e 4 }> 



In the last formula, (ip 2 ) L an d (^2)^ are obtained from (|4.53|) and given by 



(4.54) 



2 " B( " 2V F { 1)^29, 



>2v 



4(i? 2 _ 1)2 



AR{R 2 - l) 2 



F(fj) sin 29. 
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(4.55) 



(ip^) L is the stream function for the third-order Lagrangian velocity (see Appendix B): 



where D\, -D3 and D4 are constants (given by ()A.26p and ()B.4ip ). 

Thus, the Stokes drift produces a nonzero contribution to the third order outer flow given by the 
second terms on the right sides of Eqs. ()B.37p and (|B.38p . Note that the Stokes drift corrections to 
the outer flow in the lower order approximations are all zero. Also, it follows from Eqs. ()B.37p and 
(|B.38p that the Stokes drift contribution to (V"|) L vanishes in the limit R — > 00, which means that 
if there were no outer cylinder, the third order Stokes drift correction to the outer flow would be 
zero too. This can also be shown independently by treating the problem on a steady flow produced 
by an oscillating cylinder in an unbounded fluid. Thus, the appearance of a nonzero Stokes drift 
correction to the outer flow is caused by the presence of the outer cylinder. 



5. Discussion 



Equation (|4.54|) together with Eqs. (|4.4U|) . (|4,55p ~ (|4.56p represent the first two non-zero terms in 
the asymptotic expansion of the stream function for the averaged Lagrangian velocity. Let us first 
discuss the domain of applicability of formula (|4.54p . 

Domain of applicability. Our asymptotic expansion is formally valid for e <C 1 and for v = 0(1). 
In practice, we may expect that it will be valid for all e and v such that the contribution of the 
0(e 3 ) term to the right side of Eq. (|4.54p is smaller than the contribution of the 0(e 2 ) term. For 
each value of R > 1, this requirement corresponds to a domain in the space of parameters e and v. 
It is convenient to rewrite Eq. (|4.54p in the form 

^ L = e 2 [<f> (R,r) + e^$i(R,r,v,e)] sin20 + O(e 4 ) 

where <*> = sin 20 and $1 = ^ [(^) L + {V 2 ) L | f=(r _i )/e +(^) L l^=(Ji-r)/J / sin2 ^ 
Consider now the following quantity 



k(R, v , e) = t\fv max |<&i(-R, r, v, e)| / max \&Q(R,r)\, 
re[l,R] re[l,R] 

which measures the magnitude of the second nonzero term relative to the first term. We expect 
that our theory will work for all R, v, e such that >e(R, u, e) < 1, and the smaller k is, the better 
the theory should work. The level curves >c(R, u,e) = 1 in the (e, u) plane for several values of 
R are shown in Fig. [TJ For each R the domain, where we expect that the theory is valid, lies 
below the corresponding curve. Of course, the 'domains of applicability' shown in Fig. [T] are very 
approximate. Nevertheless Fig. [1] gives us an idea of where our theory might work. In particular, 
it shows that for v > 1 the interval in e within which the theory is applicable is very narrow, but it 
becomes much wider for smaller values of v. Figure [T] also shows that the 'domain of applicability' 
shrinks when we reduce the radius R of the outer cylinder. Our theory is also unapplicable to 
the case of high Reynolds numbers when the dimensionless viscosity v is comparable with e. The 
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e 

Figure 1: The curves x{R,v,e) = 1 for R — 2, 3, 5 and 10. The dashed line corresponds to v = e. For each R, the 
'domain of applicability' of the theory lies below the corresponding curve, but above the dashed line. 

dashed-line curve in Fig. [T] corresponds to v/e = 1, and the area below this curve is not covered 
by the present theory. 

Comparison with experiments. Schematic picture of the averaged Lagrangian flow is shown in Fig. 
2. The flow is symmetric relative to the axis of oscillation (the rc-axis). There are four stagnation 
points in the flow (two on the x axis and two on the y axis) and three of them are shown in Fig. 
[2j The stagnation points are symmetric relative to the x and y axes and their distance from the 
origin is denoted by r s in Fig. [2] (note that the thickness of the inner vortex structure is r s — 1). 
Examples of the streamlines produced by formula (|4.54p are shown in Fig. [3j The streamlines 
shown in Figures [3^, [3b and [3fc (which correspond to e = 0.05, 0.06 and 0.07 for v = 0.5 and 
R = 4) are similar and qualitatively consistent with the experimental observations (see, e.g., 0], 
[2l| ] ) . Figures EK, [3b and Efc show that the thickness and the intensity of the inner vortex structure 
increases with e. At some value of e the inner vortex structure becomes dominating (Fig. [3H). 
This is inconsistent with the experimental observations and corresponds to the situation where the 
0(e 3 ) term in (|4.54p is larger than (or comparable to) the 0(e 2 ) term. Our theory is not valid in 
this situation. 

Typical profiles of the azimuthal velocity are shown in Fig. [H Qualitatively the behaviour of the 
azimuthal velocity is similar to the velocity profiles measured in the experiments by Bertelsen et. 
al. 0], but there is no quantitative agreement, because in the experiments e ~ 0.3, which is outside 
the area of applicability of our theory. 

The averaged flow structure does not change for a wide range of values of R. Figure [5] shows 
the dependence of the distance of stagnation points from the origin r s on the radius of the outer 
cylinder R and on e. Qualitatively the behaviour of r s (as a function of R) agrees with experimental 
observation of Bertelsen et. al. [6J. The structure of the averaged flow changes slightly for smaller 
values of R (when R is not very different from 1). In this case, the boundary layer at the outer 
cylinder results in appearance of an additional narrow vortex system near the outer cylinder. Figure 
[6] shows the streamlines for the Lagrangian velocity for v = 0.4, e = 0.048 and R = 1.7. One can 
see a narrow boundary layer at the outer cylinder which is similar to the boundary layer at the 
inner cylinder but considerably weaker. 
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Figure 2: Typical streamlines of the averaged Lagrangian flow. 



a) e=0.05 




(b) e=0.06 




(d) e=0.08 




gure 3: The streamlines of the averaged Lagrangian flow for v = 0.5, R 
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One cylinder in the fluid that extends to infinity. 
Lagrangian velocity reduces to 



In the limit R — > oo, the stream function for the 



e 2 & + 6 3 



^3 + (^) L +0(e 



(5.1) 



where 



^2 



sin 26* , ip{ 



75 + sin 26, (^) J 



'2v 



F(Q sin 29. 



Wang |5i] studied the steady flow produced be a fixed circular cylinder placed in an oscillating flow. 
He applied the method of matched asymptotic expansions and obtained a uniformly valid expansion 
of the stream function under the same assumptions as in the present paper. His theory however 
was not directly applicable to the steady streaming flow produced by a cylinder oscillating in the 
fluid which is at rest at infinity, because his expansion of the Eulerian velocity field is not invariant 
under the appropriate change of the reference frame. In order to obtain the invariant velocity filed, 
one should consider the Lagrangian velocity which is different from the Eulerian velocity by the 
Stokes drift velocity of the fluid particles. This had been understood first by Skavlem and Tj0tta 
0], and Bertelsen et. al. jg] had corrected Wang's theory by taking account of the Stokes drift. 
However, the averaged stream function for the Eulerian velocity in Wang's theory (given by Eq. 
(3.36) in 0] and Eq. (25) in jf|) is incomplete because of the absence of the 0(e 3 ) term associated 
with the averaged outer flow. In order to compute this term, it is necessary to obtain first the 
second nonzero term in the inner expansion of the averaged stream function and then, following 
Van Dyke's recipe ^2^, match the first two terms of the outer expansion with the first two terms 
of the inner expansion. This will result in a correct 0(e 3 ) outer flow term, which then can be 
incorporated in the composite formula valid for the whole flow domain. It can be shown that if 
this is done, then the uniformly valid expansion for the averaged stream function becomes 



1 



'2v 



-75 + 



49 



sm26 + 0(e 4 



where G(£) = e ^ (e ^ + 12 cos £ + 8 sin £ + 4£ sin £) . Exactly the same formula can be obtained by 
the Vishik-Lyusternik method. If we now take into account the Stokes drift, we obtain Eq. (|5.ip . 
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This proves that (at least) up to 0(e 3 ) terms, the averaged stream function for the Lagrangian 
velocity is the same both in the oscillating (with the cylinder) reference frame and in the fixed (in 
space) reference frame. 

A remark on the results of Duck and Smith [12,]. As was mentioned earlier, the expression for ^ 
obtained here agrees with formula (3.33) in [12| (as well as with equation (3.1) in [ll|]). However, 
the boundary layer parts of the expansions are different. The most essential difference is that in our 
expansion of the Eulerian velocity there is a nonzero steady boundar y lay er at the inner cylinder 
in the first order in e, while it is present only in the second order in [12|. This is a result of the 



transformation of coordinates employed in [121 ]. The transformation maps the gap between two 
eccentric cylinders onto the annulus between two cylinders whose axes coincide and it is a time- 
dependent transformation. In order to obtain our expansion, one needs to do the following: (i) to 
write down a composite expansion using the formulae for inner and outer expansions obtained in 
[l2l | . (ii) to perform the inverse transformation of the coordinates, and (iii) to expand the resulting 
velocity field in Taylor's series in e. The result will almost certainly coincide with our expansion 
of the averaged stream function for the Eulerian velocity up to 0(e 2 ) terms. 

Remarks of the steady streaming theories for high Reynolds numbers. There are quite a few papers 
dealing with the steady streaming in an unbounded fluid produced by an oscillating circular cylinder 
(see 0, 0, @, Qj 0] ) a t high Reynolds numbers such that Re s = 0(1) or Re s S> 1 (in our notation 
this means that v/e 2 = O(l) or v/e 2 -C 1, respectively). In all these papers, coordinate systems 
oscillating with the cylinder are used, i.e. what is being really solved is the problem about a steady 
streaming flow produced by the fixed cylinder placed in an unbounded oscillating flow. As the 
present study shows, the problem with a fixed cylinder is not equivalent to the problem with an 
oscillating cylinder not only within the boundary layer but also in the outer flow. In fact, the 
situation is even more difficult, because for v/e 2 = 0(1) (and even more so for u/e 2 <C 1) it is 
unclear how to solve the problem in the frame of reference fixed in space. In the present study, 
we employed the standard trick: using the fact that the amplitude of oscillations was small we 
expanded the solution in Taylor's series about the averaged position of the oscillating cylinder and 
thus transferred the no-slip boundary conditions at the moving surface of the oscillating cylinder 
to the fixed surface of the cylinder at its averaged position. This allowed us to formulate our 
asymptotic expansion in terms of a sequence of boundary value problems with fixed boundaries. 
The thickness of the boundary layer on the inner cylinder that appeared in our expansion is 0(e), 
which is much larger than the 0(e 2 ) displacement of the inner cylinder from its averaged position, 
and this is what justifies the transfer of boundary conditions from the moving boundary to the 
fixed one. If however we considered the case of v = e 2 (Re s = 1), the thickness of the Stokes layer 
would be 0(e 2 ), i.e. of the same order as the displacement of the cylinder, and therefore it would 
be impossible to justify the transfer of boundary conditions. Even if we worked in the frame of 
reference fixed with the oscillating cylinder, we would encounter the problem of calculating the 
Stokes drift velocity in the boundary layer because of non-analytic dependence of the boundary 
layer velocity on e. 



In ll| and [12J, the same problem as in the present paper, i.e. the steady flow between two 
cylinders produced by small-amplitude oscillations of the inner cylinder, had been studied in the 



case of Re s = O(l) (and Re s 3> 1). Haddon and Riley [ll| had realised the impossibility of the 
transfer of the boundary conditions from the moving surface to a fixed one in this flow regime and 
addressed the problem by employing two different coordinate systems for boundary layers on the 
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Figure 6: The streamlines of the averaged flow for v = 0.4, e — 0.048 and R = 1.7 



inner and outer cylinders. In the end, however, they had applied the transfer of the boundary 
conditions at the inner cylinder for the steady outer flow, and it is unclear whether this can be 
justified. In [12], a conformal mapping that maps the gap between two eccentric cylinders onto 
the annulus between two cylinders with a common axis was employed. The subsequent analysis 
had been done using the transformed coordinates. Neither the inverse transformation to physical 
coordinates, nor the Stokes drift had been computed, and as the above discussion indicates, these 
are the questions where potential problems may arise. Thus, in spite of a considerable progress in 



this area and a good agreement with experimental results achieved by the theory (see, e.g., 121]). 
there are still certain unanswered questions concerning the steady streaming at high Reynolds 
numbers. 
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6. Appendix A. Explicit expressions for the right hand sides of Eqs. (13.181) (13. 20|) . 
(|jOID - (!3H3D , (OOH) and (|OS|) 

Functions F fc a , G a k and H% for fe = 1,2, 3 in Eqs. (^18l) - (l3T20l ) are given by 

Fi = -«jlr=i«oe + ^Se + ^ (A-l) 

n = -4U<e-KU+e4X=i)^-^(^+4l,=i<) 

- 4\ r=1 v a + v « c - + v a m - v%) + - fp a oe , (A.2) 
-de (vl\ r=1 vt + {v\\ r=1 Hvl r \ r=1 ) v a + v a Q vt - - U\ r =i wo) 

-<VU - ufvfe - 4rlr=l < ~ «olr=l V l ~ ( U Wr=l +f«0rlr=l) V ~ V o\r=l < 

~<v a - e4l,=i vo + - + fPoe 
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+v(v% - + fv^ + v a we - 2iv a m - vf + + 2u a 0l 



oe 1 1 



(A.3) 



G? 


= 0, 




(A.4) 


Gf 






(A.5) 




= - Kt - ^lee) - u olr=l «0£ " n o 


9 \ r=1 «g + «) 2 + 2«<| r=1 wg + v (ufc - 2zfo) 


, (A.6) 








(A.7) 


^2 






(A.8) 


#3° 






(A.9) 



Functions F b , G h k and H b for k = 1, 2, 3 in Eqs. ([3T2Tjl - (|3T23l ) are given by 
P fe "..6 »7 *• 



F b = 



R v 0v R2 Poe 

'"3 = n olr=iJ u 2?y + { u l\r=R ~V u 0r\r=R) 4t] + ^ M 2lr=_R ~~ r ? n lrly=fl +"^" n 0rrlr=_R^ u 0»y 

1 " ,j| r=fl ^ b + (v\\ r=R - V vl\ r=R ) v b + v b v b + ± + v\ 

r=R 4 + {4\ r=R -v4r\r= R ) ^0 + «olr=Jl 4 1 '" 
,2 



-^0 n', 
1 

+^ + 4v b 0v - 



-i-T) 6 - ^-n b 

R2 Pie R3 Poe, 



2 ■+t;5| r=Ji ug 



G\ = 0, 
Co = 



\\ r=R u b + u b v b 

2 3 

V b V b 



u 0t uu 0rm ' 



G\ = 



u\ T — vu\ — u l Q 



\r]ri ~ u o\r=R n 0»? + u 0e\r=R V 





4 ((<*)' 


+ H\r= 


R v b o 


H\ = 


R 2 v oe 






H b = 


R 2 v w 


R 3 v oe 




Hi = 


V v b 


^ v b 





f -.6 



V 



~ R U2 ~ m Ul ~ R? U °- 
Functions Q b k and S\, in Eqs. (|3.25p and (|3.26p are given by 

Qi = 0, 
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(A.10) 



(A.ll) 



(A.12) 

(A.13) 
(A.14) 



(A.15) 
(A.16) 
(A.17) 
(A.18) 

(A.19) 



Q 2 = -cos0/(^ r | r=1 +<^ =o J + sin0/ufo| r=1 (A.20) 
Q 3 = - cos 6> / (V lr ,| r=1 +«i^| ?=0 ) + sin / (V ie | r=1 +ng e | ?=0 

-sin0/ (4| r=1 +<| ?=0 ) - cos 2 #y<^ =0 -sin0 cos 0/ 2 < ? | €=o , (A.21) 

51 = -cos0f(r)v^ =Q , (A.22) 

5 2 = - cos 6 f (v l 0r \ r =i + v i^=o) + sin 9 f ( u oelr=i +*4ls=o, 

+ sin0/u* o | r=1 -cos 2 ^^| €=0 , (A.23) 



S 3 = - cos / (<| r=1 +f2 € | ?=0 J + sin / (wiflUi -H&| €= 
+ sin^/ («l| r=1 +ng| €=0 ) - cos 2 ^^| 5=0 

+ sin cos 9 f 2 v^ =0 - sin 2 9 f —v a Q ^ =Q - cos 3 0^<^ =o . (A.24) 

Explicit formulae for constants C\, C 2 , C3 and C4 in Eq. (j4.40j) are 

Ci = R 6 (R 2 + 2) + 2R 2 + 1, C 2 = -R 8 - R 2 , 

C 3 = -R 4 (2R 2 + 1) - R 2 - 2, C 4 = i? 4 + 1. (A.25) 

Constants Z?i, D 2 , D 3 and -D4 in Eq. (|4.48j) are given by 

D 1 = -(75R W -75R 9 + 138R 8 -158R 7 

+ 106.R 6 - 134i? 5 + 131.R 4 - 79R 3 + 80R 2 - 40R + 40), 
D 2 = R 2 (53R 8 -53R 7 + 32R 6 -MR 5 + 39R A - 51R 3 + 52R 2 — 14i2 + 14), 
D 3 = 2(75i? 8 -77R 7 + 58i? 6 -7AR 5 + 59R 4 - 52R 3 + 53R 2 -40i? + 40), 
£> 4 = -(75i? 6 -79i? 5 + 31 J R 4 -43i? 3 + 44# 2 - 40i? + 40) . (A.26) 

7. Appendix B. Calculation of the averaged Lagrangian velocity 

In polar coordinates, Eq. (|4.49[) is equivalent to 

r T = e 2 [u(r,9,T)-r t ]. (B.l) 
r9 T = e 2 [v{r,6,T) - r0 t ] . (B.2) 

We already know that 

u = 4 + e (u\ + ul + ul) + e 2 (u 2 + < + uj) + • • • , (B.3) 
« = «j + + + e («! + u? + uj) + e 2 (t£ + + (B.4) 

We seek the solution of (jB.ip and (|B.2|) in the form 

r = r + en + e 2 r 2 + . . . , 6* = + e0i + e 2 6» 2 + . . . (B.5) 
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On substituting ()B.3[) ()B .5|) in (jB.ip and (|B.2p and collecting terms of equal powers in e, we obtain 
the following sequence of equations: 



and 



ror = 0, r <9 T = 0, 

rir = 0, r 9 lT + ri9 0T = 0, 



d T r n + d t r n _ 2 = U n -i, 

n n—2 

^ r m d T + ^2 r m dt9 

n—m—2 = Vn-2, 



(B.6) 
(B.7) 



(B.8) 
(B.9) 



m=0 



m=0 



for n = 2, 3, . . . Functions J7 n and V n are obtained by substitution of (|B.3[) — (|B.5j) in the right sides 
of Eqs. (|B.ip and (|B.2j) and subsequent Taylor's expansion of all terms. There is a subtle technical 
trick here. Functions which depend on boundary layer variables £ = (r — l)/e and rj = (R — r)/e 
are not analytic at e = 0. Therefore, for these functions we use Taylor's expansions at r = ro + er% 
rather than at r = ro + er\. For example, 

ug(£,0,T) = ^^^-+ri + er 2 + ...,0o + e^i + ---,r 

= <(£* ,9 , t) + e [r a 4;(£*, O , r) + 0i ufo(f A, r)] + . . . 
r\. Below we use the following notation: (•) | = (•) \ T=TQ g= g and 

or (■)!.= (■) 



where £* 



(') I* (') I ^ = ro_A +r lj o=e 



depending on whether the quantity being evaluated is a function of £ or rj. The explicit formulae 
for U„ and V„ are 



Vi 

lh 



Vo 



"olo. y o = ^olo + v ol* +^ol*> 
[t4+rit4r + ^i«otf] lo +«ol* +«oL> 

[vj + n^ r + Oivio] | + [u? + r 2 v^ + 0^4] |, + [v\ - r 2 v b 0v + 9 lV b oe 
u\ + riui r + 9iu[ + - (r?u 0rr + 2r 1 9 1 u i 0rd + 0?u O 0<?) + r 2 u 0r + 6» 2 u oe 



+ [«? + r 2 u^ + ffitig, 



06»J 



+ 



u\ - r 2 u b 0v + 6»iUo e 



(B.10) 
(B.ll) 

(B.12) 



(B.13) 



v\ + nv[ r + 9 x v\ e + - (r 2 x vi rr + 2n0it;5 rfl + 0??4 e ) + r 2 «j r + 



+ 



+ 



t> 2 + r 2 i;% + 9 lV %- {r 2 2 v a m + 2r 2 9 lV ^ e + 9 2 v% 69 ) + r 3 vfc + 6 2 vfo 



1 



4 ~ r 2 v° lr , + 9 lV ° ie + - I rf^ - 2r 2 9 lV ° 0rie + 



q2„.6 



+ 02Vm 



(B.14) 



It follows from Eqs. ()B.6p and (|B.7p that ro, 9$, r\ and 9\ do not depend on fast time r, i.e. 
r o = r o(^)i 00 = 0o(i), r i = r i(t) and 0i = 9\{t). Equations (jB.lOp and the fact that u = 0, 
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i>5 = U, Vq = 0, Vq = imply that Uq = and Vq = 0. Then, averaging Eqs. (|B.8[) and (jB.9[) 
(for n = 2), we find that rot = and ro#ot = 0. Therefore, ro = const and 0o = const, and we 
use ro and 0o to identify fluid particles: (r(ro, 9q, t, r), 0{tq, 0q, t, r)) represents the trajectory of a 
fluid particle whose averaged (in r) position at t = was (ro,#o)- From Eq. (jB.lip . we deduce 
that U\ = 0. Averaging Eq. (|B.8[) for n = 3, we find that r\t = 0, so that n = const. We choose 
T\ = 0. Note that this implies that (•) |*= (•) | in Eqs. p^ - ([R14]) . Averaging Eq. (|R9]> for 
n = 3 yields 

(B.15) 



«1 + r 2 < 



0€ 



To proceed further, we need to compute r2V^ and T2Vq v 
part of T2 does not make any contribution to these two terms, i.e. r2V^ 
r^Uojj • The oscillatory part of T2 satisfies the equation 



. First we observe that the averaged 



r 2V b o v 



>'2t 



U 



it. 



olo= Uo{r ,0 ,T) 



(which is simply the oscillatory part of Eq. (|B.8j) for n = 2). Substitution of (|4.ip yields 



r 2 



i? 2 - 1 



i? 2 



1 cos 6q. 



(B.16) 



This implies that 



r 2 v^ 



cos 6<o / i? 



i? 2 -l 



COS fo 
i? 2 -l 



i? 2 



0/;- 



On the right sides of these equations we have functions and Vq^ which are nonzero only within 
boundary layers with thickness 0(e). It is therefore natural (and consistent with our asymptotic 
expansion for the velocity) to replace ro by 1 + e£o in the first equation and by R — erjo in the 
second equation (where £o = ( r o — l)/ e an< ^ % = (R ~ r o)/ e ) an d expand everything in Taylor's 
series. This leads to the expressions 



r 2 v. 



0£ 



1 



2£oR 2 . m 2- 

e—n h 0(e 

R? - 1 v ' 



COS 0Q f 

e R{&-l) cosg o7< + Q( g 2 ) 



(B.17) 
(B.18) 



Now we substitute 0(1) terms from (|BT7]) and (fB~18l) in Eq. (fB~Tf)j) and move higher order terms 
to appropriate higher order equations. Equation (|B.15p becomes 



roVit 



Vf + COSQ / Vq; 



0£ 



The first term on the right side of this equation represents the averaged Eulerian velocity, while 
the second term is associated with the Stokes drift. Finally, substitution of (|4.15p reduces this 
equation to ro0\t = 0, which, in turn, implies that 0\ = const. Again, we choose 0\ = 0. Thus, 
the averaged Lagrangian velocity has no 0(e) boundary layer term. Averaging the equation for 
yields 

" ' 1 f 1 1 " (B.19) 



T2t 



h + r 2 u 0r + V2 U 06 



+ 
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-r 2 u b 0v 



Now we need an explicit formula for 62- The equation for 82 is 



12t = — \ V \ +V Q \ Q +V Q \ Q 

On substituting Eqs. (|4.ip . (|4.4p . (|4.8p and integrating over r, we obtain 

9~2 = 62 + 82+ ®2 

Si sin0 o 1 (R 2 \ 



2R 2 sinO, 
R 2 -l 
2 sin 6q 



R(R? 



-Re[Ce-^ 0+iT j + 0(e), 
-Re (Ce"™ +iT ) +0(e). 



Now we return to Eq. (|B"T9]) . It follows from (gJJ), g^D, (jOjl . (jB~T6|) . (|R2H - (|R231) that 



r 2 u l 0r 



7 2 "06» 



0, e\u i oe 



o, 2 b <, 



O(e), f 2 u b 0ri =0(e), 



sm9 fv% + O(e), f 2 ^ q = -cos0 o /^ + O(e). 



After substitution of (1R24]) and ([R251) in Eq. (fB~19j) . it reduces to 

= uKro,^). 



(B.20) 
(B.21) 

(B.22) 
(B.23) 

(B.24) 
(B.25) 

(B.26) 



Note that the right side of this equation is different from what we would have if we used the 
averaged Eulerian velocity: the second order term in the expansion of u is u\ + u\ rather than u 2 - 
This means that the Stokes drift kills the boundary layer term. 

Now let us deduce equation for 62- On averaging the equation for #4 and using (|B.14p . we obtain 

1- 



v l 2 + r 2 v l 0r + 6 2 v l QB 



+ 



v 2 + r 2 v^ + 2 r 2 u 0££ + r 3<£ + d ^ v k 



+ 



1- 



v\ - r 2 v\ n + ^r%v b 0vv - r 3 v b 0j] + 8 2 v b v 



r 2 8 2r + T 2 , 



where all quantities on the right side of this equation are evaluated at (ro,#o) an d where 

2% 



COS UQ-p^ -J V% f + COS 8q- 



represent the contribution that comes from O(e) terms in (IB.17P and ()B.18p . To proceed further, 
we need an explicit formula for r 3 . We have the equation r 3r = u\ | +Uq | +Uq | , from which, on 
substituting (|4.6p . (|4.9|) . (|4.12p and integrating over r, we obtain 



(B.27) 



(B.28) 



rn + rq + r 



3' 



2R 



(R 2 - l) 2 
2i? 2 



R + l 



C 



R 3 + l 



Be [ —e tT ) cos f) (] . 



R(R 2 - 1 



■i.'r I - e 

7 



-7770 -Ht 



COS &Q. 



(B.29) 
(B.30) 

(B.31) 

(B.32) 
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After tedious but elementary calculations with the help of (jRPol) . (lR20l) - (lR23l) . (fB39j)HfR32j) . 
Eq. (|B.27j) can be reduced to 

roht = 4 + ^2+^2 + ^2 
where v 2 represents the Stokes drift velocity of the fluid particles and is given by 



v 2 



cos Vq 

1 



fv^-fv^ 



+ 



2 cos 6>o 
R 2 -l 



R 2 



9 v o$ ~ j ^ 



C 



1 



+ 



R 



Vo f v o v - Re 



-fijo+ir 



(B.33) 



Here g(r) = Re(Ce lT /j). In order to obtain the 0(e 3 ) term in the expansion for the stream 
function, we need to compute the components of the averaged third-order Lagrangian velocity in 
the outer flow. Averaging the equations for r§ and #5 and ignoring all boundary layer terms, we 
obtain 



ro4 



u 3 + r 2 U lr + r 3 U 0r + °2 U W + °l U 



00' 



-r^3r ~ r 3»2r + ^3 + % + r 3 V 0r + ^16 + 



(B.34) 
(B.35) 



All terms on the right sides of these equations are evaluated at {tq,9q). Now we need to find 9\. 
The equation for 9\ is obtained from (IB.9P for n = 3 by taking its oscillatory part and ignoring all 
boundary layer terms: 

r 9i T = v\(r ,e ,T). 
Inserting (4.17) into this equation and integrating over r, we find that 

» 2i? 



(R 2 - l) 5 



R 3 + I s 

(R + l) + =^—)g(T) sm9. 



(B.36) 



Substitution of (|Q]L gH|), (Q5j| . (jRl6]) . (jB~2Tj) . (|BT30]l and (jB~36|) in Eqs. (jB~34|) and 

(lR35l) yields 



4i? 



r 3< 



^3 + 



(i? 2 - l) 2 

AR 

(r 2 - iy 



f'g 



f'g 



cos 29 



sin 29 



Further calculations yield 



'2v 



AR{R-l){R 2 - l) 4 
V2^ 



where 
Hence, 



AR(R-l)(R 2 - l) 4 

D 2 = D 2 - AR 2 (R - 1)(R 2 - l) 2 . 

t>2 



D 1 D 2 A 
— + ^- + D 3 r + Dy cos 29, 
r r3 J 

^-D 3 r- 2D 4 r 3 | sin20, 



2^ 



8R(R-1)(R 2 - l) 4 



£>l + ^| + A^ 2 + sin 20 



(B.37) 
(B.38) 

(B.39) 
(B.40) 

(B.41) 
(B.42) 
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